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Abstract. - In recent years there have been a number of proposals to utilize the specificity of 
DNA based interactions for potential applications in nanoscience. One interesting direction is 
the self-assembly of micro- and nanoparticle clusters using DNA scaffolds. In this letter we 
consider a DNA scaffold method to self-assemble clusters of "colored" particles. Stable clusters 
of identical microspheres have recently been produced by an entirely different method. Our DNA 
based approach self-assembles clusters with additional degrees of freedom associated with particle 
permutation. We demonstrate that in the non-equilibrium regime of irreversible binding the 
self-assembly process is experimentally feasible. These color degrees of freedom may allow for 
more diverse intercluster interactions essential for hierarchical self-assembly of larger structures. 



DNA has attracted significant attention for its poten- 
tial applications in nanoscience ( [1], [2], [3], [4], [5], [6], 
[7], [8], [9], [10]). One recent non-DNA based advance 
is the self-assembly of stable clusters composed of iden- 
tical microspheres [11]. In this letter we consider the 
self-assembly of micro- and nanoparticle clusters similar 
to those of [11], where DNA scaffolds govern the self- 
assembly process. The plan for the letter is the following. 
We first introduce the basic strategy of our self-assembly 
proposal. The goal is to maximize the yield for a particu- 
lar type of cluster we call the star cluster. We analytically 
calculate the yield of the star cluster in the regime of ir- 
reversible binding. The analytical results are compared 
to the numerical results for the full aggregation equations. 
From an experimental perspective, the most important re- 
sult is the determination of an optimal concentration ratio 
for experiments (see Eq. 8). To conclude we discuss the 
experimental feasibility of the self-assembly proposal. 

The basic idea behind the procedure is as follows (see 
Fig. 1). Particles are functionalized with single-stranded 
DNA (ssDNA) markers which determine the particle color. 
There may be many DNA attached to each particle, but on 
any given particle the marker sequence is identical. One 
then introduces DNA scaffolds to the system. The scaffold 
is a structure with / ssDNA markers, each marker com- 
plementary to one of the particle colors. Hybridization of 
the ssDNA markers on the particles to those on the scaf- 



fold results in the formation of colored particle clusters. 
Because there are many DNA attached to each particle, 
clusters can form which contain more than one scaffold. 
The essential goal of the procedure is to maximize the 
concentration of a particular type of cluster which we de- 
note the star cluster. The star cluster contains one and 
only one scaffold to which / particles are attached, each 
particle having a distinct color. 

We should note that the role of the scaffold could also be 
played by a patchy particle ( [12], [13], [14]). For example, 
these patches are regions on the particle surface where 
one can graft ssDNA markers. In this case there may 
be several DNA connections between a patch and colored 
particle. Our conclusions will still be valid, provided the 
patch size is chosen so that a patch interacts with at most 
one particle. 

Previously we performed an equilibrium calculation to 
determine the yield of the star cluster [15]. The results 
of that study indicated that the concentration of scaffolds 
must be kept very small to prevent the aggregation of 
larger clusters. From an experimental perspective this 
result is somewhat disappointing, since the overall yield 
of the star cluster is proportional to the scaffold concen- 
tration. The situation is considerably improved in the 
regime of irreversible binding of particles to scaffolds. In 
what follows we present a calculation for the yield of the 
star cluster far from equilibrium. 
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the star mass fraction M* = (/C*)/Ct, 



DNA scaffold 



Colored particle 



Fig. 1: DNA scaffolding. A graphical depiction of the 
scheme for self-assembling star clusters using DNA scaffolds. 
In the diagram (not drawn to scale) the scaffold funcionality 
/ = 4. 



To understand the basic physics behind the aggregation 
process we consider the mobility mismatch between the 
particles and the scaffolds. In solution, a particle with 
radius R ~ lpm has a diffusion coefficient given by the 
Stokes-Einstein relation D — ksT /6ttt]R. On the other 
hand, the size of the scaffold a ~ lOnm. As a result 
the scaffolds diffuse R/a ~ 100 times faster than the par- 
ticles. To first approximation the resulting aggregation 
is a two stage process. In the first stage the particles 
recruit different numbers of scaffolds via the fast scaffold 
diffusion and subsequent DNA hybridization. Since we 
consider the regime of strong binding where these bonds 
are irreversible, the result is a Poisson distribution over the 
concentration of particles with m scaffolds attached. Let 
Ci denote the concentration of particles with color i, and 
c denote the total concentration of scaffolds. The total 
particle concentration Gtot = Y^J—x Ci- The concentra- 
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of particles of color i with m scaffolds attached 
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In the second stage there are no free scaffolds left in 
solution, and these particles decorated with scaffolds ag- 
gregate to form the final clusters. The seed to build a star 
cluster is a particle of any color with exactly one scaffold 
attached. This seed must aggregate with / — 1 particles 
of different colors, each of which has no scaffolds. We now 
calculate the concentration of the star cluster C*. The 
yield of the desired star cluster is quantified in terms of 



M* = 
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Here x = fp is the scaffold functionality / multiplied by 
the concentration ratio p. By choosing p = 1/ f the mass 
fraction attains a maximum of exp(— 1) ~ 0.37! This re- 
sult indicates that by selecting the appropriate scaffold 
concentration, in the noncquilibrium regime up to 37% of 
the particles will aggregate to form star clusters. This is 
a significant improvement over the situation in the equi- 
librium regime. 

This treatment of the problem captures the physics of 
star cluster formation, but it does not account for the loss 
of star clusters due to aggregation. In particular, as long 
as there are scaffolds with markers available for hybridiza- 
tion, when these scaffolds encounter a star cluster they 
can aggregate to form a larger cluster. We now estimate 
how this aggregation effects the final concentration of star 
clusters. 

Consider the beginning of the second stage in our aggre- 
gation process. There are no longer any free scaffolds in 
solution, but a scaffold can have up to / — 1 DNA markers 
still available for hybridization. We would like to deter- 
mine how the star cluster mass fraction M* (y) changes as 
a function of the fraction of saturated scaffolds y. Here 
a saturated scaffold has particles hybridized to all / of its 
DNA markers, and is therefore unreactive. If s is the 
expectation that a slot on the scaffold is filled, then the 
fraction of saturated scaffolds is y = s* _1 . The average 
number of open slots on a scaffold is (/ — 1)(1 — s). Con- 
sider filling an open slot on the scaffold. The probability 
that the particle which filled the slot was part of a star 
cluster is M*{y). The average rate r(y) at which star 
clusters are lost to aggregation is then 

r(y) = ~MM^ i(f - - *)] = M*{y)y- a . (4) 

Here the exponent a — (/— 2)/(/ — 1). We can then con- 
struct a differential equation for M„ taking into account 
this loss due to aggregation. 
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In the absence of this loss term the result of the calculation 
should recover our previous result Eq. 3. This zeroth 
order approximation is just Mi°\y) = xy exp(— xy) which 
gives the correct star cluster concentration once all of the 
scaffolds are saturated (y = 1). To simplify the analysis 
a bit we take a = 1 which is an excellent approximation 
in the limit of large scaffold functionality /. This is an 
inhomogcncous first order differential equation which can 
be solved by introducing an integrating factor u{y) = y x . 
The initial condition which must be satisfied is M*(0) = 0. 
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We are interested in the final star mass fraction M* , which 
isM*(y = l). The result is 



fc=0 



(-*) 
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fc + 1 x + k + 2 
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= sexp(-a;) + x 2 E^ x {x) - x 1-i T(l + x) 
Here r(x) is the gamma function and E v (x) = 

oo 

J t~ v cxp(—xt)dt is the exponential integral of order v. 
1 

We can perform a similar type of analysis in the case 
when there is only one particle color. In this case the / 
ssDNA markers on the scaffold all have identical sequences 
complementary to this color. It turns out that the result 
for the mass fraction is the same. Because the mass frac- 
tion is the same in both cases, we can gain insight into the 
behavior of the system with many colors by analyzing the 
much simpler one color system. To test our predictions, 
we numerically solved a system of differential equations 
which models the irreversible aggregation between parti- 
cles (one color) and scaffolds. 
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Fig. 2: Star cluster mass fraction. The mass fraction M* 
as a function of p for scaffolds with functionality / = 3 (red) , 
4 (blue), 5 (green), and 6 (black). The results determined 
numerically from the full solution of the Smoluchowski coag- 
ulation equation (markers) can be compared to results of the 
anlaytical calculation (lines) Eq. 6. The resulting agreement 
is good up to an overall normalization factor 7/ in the range 
Cij S ' KjjjjCjj (7) 1.2 to 1.5 which normalizes the analytical curves. 
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This equation is the Smoluchowski coagulation equation 
[16] adapted to our system. Cij is the concentration of 
the cluster with i scaffolds and j particles, ifyj'j' is the 



rate constant for the irreversible reaction CV 



a, 



i'j' 



Ci+iij+j> . We assume that the rates are diffusion limited 
in which case we can estimate the rate for any pair of 
clusters by Kiji>ji = 4irD s Ri. The larger cluster with 

1/3 

hydrodynamic radius Ri ~ n ; plays the role of a sink. 
Here n\ is the number of particles in the larger cluster 
and D s = ksT '/QnriR s is the diffusion constant for the 
smaller cluster. To simply matters we only consider tree 
like structures, i.e. we do not consider the formation of 
clusters with internal loops. We have truncated the set of 
equations by considering clusters with a maximum of 10 
scaffolds. 

By solving these equations we can determine the con- 
centration of stars C* = C\f in this notation and test the 
validity of our two stage ansatz. As indicated in Fig. 2, 
the result of our analytical calculation matches the results 
of the full numerical calculation up to an overall normal- 
ization factor of order unity. Several points are in order. 

The optimal concentration ratio p for experiments is 
easily determined from — 0. The result is £ m ax — 
0.47. For scaffolds of functionality / the concentration 
ratio should be chosen as: 



P 



0.47 



(8) 



Note that the maximum attainable star cluster yield 
M*(x max ) — 1/4 does not decrease with increasing /. In 



fact, the numerical results predict a slight increase in star 
cluster yield for larger /. Solving the aggregation equa- 
tions becomes computationally expensive, but it can still 
be done by reducing the maximum number of scaffolds in 
a cluster. For example, considering clusters with up to 
5 scaffolds for / = 10 gives M*(x nmx ) ~ 0.3. These re- 
sults are important from the perspective of experimental 
feasibility for the self-assembly method. This is to be 
contrasted with the earlier equilibrium treatment. There 
the condition to suppress the aggregation of larger clus- 
ters imposed a fairly strict constraint [15] on the concen- 

tration ratio p < J 1 / 2 (j) . From the perspective of 
self-assembling stars with large / this renders the regime 
of irreversible binding far more appealing than the equi- 
librium regime. 

If an experiment is performed with the optimal concen- 
tration ratio, the clusters which self-assemble are easily 
separated by density gradient centrifugation [17]. In this 
regime most of the particles are monomers, in star clus- 
ters, or in saturated two scaffold clusters. These clusters 
contain, 1, /, and 2f — 1 particles respectively. The dis- 
parity in hydrodynamic radius and sedimentation velocity 
of these clusters makes the separation procedure experi- 
mentally feasible. 

In this letter we considered a DNA scaffold method for 
self-assembling star clusters of / colored particles. By 
taking advantage of the mobility mismatch between par- 
ticles and scaffolds, we were able to formulate a nonequi- 
librium calculation of the star mass fraction. The re- 
sults of the calculation were compared to the numerical 
results of the full Smoluchowski coagulation equation for 
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the system. Good agreement is established between the 
analytical calculation and the numerics. In the regime of 
irreversible binding the yield of the desired star cluster is 
drastically improved in comparison to earlier equilibrium 
estimates. In noncquilibrium we find an experimentally 
feasible regime for the self-assembly of star clusters with a 
maximum mass fraction ~ 1/4. We determined the opti- 
mal concentration ratio for an experimental implementa- 
tion of our proposal. The additional color degrees of free- 
dom associated with particle permutation in these clusters 
makes them ideal candidates as building blocks in a fu- 
ture hierarchical self-assembly scheme. In addition, these 
clusters can serve as the starting point to self-assemble 
structures of arbitrary geometry [18]. The experimen- 
tal realization of self- assembling star clusters using DNA 
scaffolds would constitute an important step towards re- 
alizing the full potential of DNA mediated interactions in 
nanoscience. 

This work was supported by the ACS Petroleum Re- 
search Fund (Grant PRF No. 44181-AC10). We acknowl- 
edge L. Sander, G. Ghoshal, and B. Karrer for valuable 
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